Observation of discrete gap solitons in binary waveguide arrays 
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We report on the first experimental observation of discrete gap solitons in 
binary arrays of optical waveguides. We observe the soliton generation when 
the inclination angle of an input beam is slightly above the Bragg angle, and 
show that the propagation direction of the emerging gap soliton depends on 
the input power as a result of an inter-band momentum exchange. 
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Nonlinear periodic photonic structures such as arrays 
of optical waveguides attracted a lot of interest due to 
the unique ways they offer for controlling light. Periodic 
modulation of the refractive index breaks the transla- 
tional invariance and produces effective discreteness in 
a continuous system, opening up many novel possibili- 
ties for manipulating light propagation, including light 
localization in the form of discrete optical solitons 1 ' 2 . 

Recently, it was suggested that a novel type of discrete 
optical solitons, the so-called discrete gap solitons, can 
be generated in the binary waveguide arrays, by a single 
inclined beam 3 or two input beams 4 . The binary arrays 
of optical waveguides are specially engineered photonic 
structures consisting of periodically alternating wide and 
narrow waveguides [see an example in Fig. 1(a)]. Dis- 
crete gap solitons in such structures can be considered 
as a nontrivial generalization of the spatial gap solitons 
recently observed in the waveguide arrays 5 and optically- 
induced photonic lattices 6 . Gap solitons in binary ar- 
rays are associated with the fundamental modes strongly 
confined in narrow waveguides, which lowers the power 
threshold for the soliton excitation, and may result in 
reduced radiation losses compared to the gap solitons 
based on radiation or higher-order bound modes in the 
conventional arrays. In this Letter, we report on the 
first experimental observation of discrete gap solitons in 
binary arrays of optical waveguides and demonstrate a 



novel method of efficient steering of gap solitons based 
on the inter-band momentum exchange. 

We investigate spatial beam self-action and soliton 
formation in the fabricated etched arrays of the 5 mm 
long AlGaAs waveguides with the effective refractive in- 
dex contrast 0.0035 (see Ref. 2). The binary arrays 
are made of wide (4 jam) and narrow (2.5 jam) waveg- 
uides with 4 /am edge-to-edge spacing, and accordingly 
the full period is d =14.5 jam as illustrated in Fig. 1(a). 
The beam propagation in this structure can be de- 
scribed by the normalized nonlinear Schrodinger equa- 



tion, i 



OE 

dz 



• v(x)E + \E\ 2 E = 0, where E(x, z) is 



the normalized envelope of the electric field, x and z are 
the transverse and propagation coordinates normalized 
to the characteristic values x s = 1/am and z s = 1mm, 
respectively, D = z s X/ (Annox 2 ) is the beam diffrac- 
tion coefficient, no = 3.3947 is the average medium re- 
fractive index, A = 1.5 jam is the vacuum wavelength, 
v(x) = 2An(x)7rno/A, and An(x) is the effective modu- 
lation of the optical refractive index. 

In the waveguide arrays, the wave spectrum possesses 
a band-gap structure with gaps separating bands of the 
continuous spectrum associated with the spatially ex- 
tended Bloch waves of the form, E(x,z) = ip(x) ex.p(ikx+ 
i/3z), where the propagation constant /3 is proportional 
to the wave- vector component along the waveguides. In 
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Fig. 1. (a) Normalized refractive index in a binary waveg- 
uide array (gray) and Bloch-wave profiles at the top of 
first (solid) and second (dashed) bands. (b,c) Disper- 
sion of Bloch waves and the Bloch-wave excitation coeffi- 
cients vs. the transverse wave- number, for the first (solid) 
and second (dashed) bands. Fourier spectrum of an in- 
put beams with 20% larger then Bragg incident angle is 
shown in (c) with shading, (d) Characteristic profile of 
a discrete gap soliton in the binary array. 



Fig. 1(b), we plot the Bloch-wave dispersion curves and 
mark two types of spectrum gaps appearing due to (i) to- 
tal internal reflection and (ii) Bragg-type scattering. 

When the medium nonlinearity is self-focusing, the dis- 
crete gap solitons may appear near the upper edge of 
the second band 4 . One of the distinguishing features of 
such solitons is their localization at narrow waveguides, 
as shown in Fig. 1(d). As was suggested in Ref. 3, the 
gap solitons can be excited by an input Gaussian beam 
which is inclined slightly above the Bragg angle. Each 
Fourier component of an incident beam excites a super- 
position of Bloch waves belonging to different bands, ac- 
cording to the values of excitation coefficients C n plotted 
in Fig. 1 (see details in Ref. 7). If the input beam spans 
several periods of the underlying lattice (we choose the 
45 jam FWHM beam) and its spectrum is narrow, then 
the excitation of the second band becomes dominant [see 
Fig. 1(c)], facilitating the formation of a gap soliton. Al- 
though fully controlled generation of gap solitons can be 
achieved only in a two-beam excitation scheme 4 6,8 , the 
single-beam approach is simpler to implement. Addition- 
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Fig. 2. Experimentally measured output intensity distri- 
bution vs. the inclination angle of the input Gaussian 
beam. 



ally, the simultaneously excited Bloch waves of the first 
band move in the opposite direction, resulting in "self- 
clearing" of the gap soliton, and allowing to control the 
soliton velocity by changing the input beam intensity 3 . 
We note that all Bloch-wave components tend to trap to- 
gether into a multi-band breather if a normally incident 
narrow beam is used 9 , and complex Bloch-wave interac- 
tion occurs for the inclined narrow beams 7 . 

In our experiments, first we study the linear beam 
propagation (see Ref. 2 for description of experimental 
set-up). The dependence of the output field distribution 
on the input angle is shown in Fig. 2. We note that, in a 
binary waveguide array, the Bloch waves corresponding 
to the first and second bands can be easily distinguished, 
since they are primarily localized at wide and narrow 
waveguides, respectively. We find that the band-2 exci- 
tation becomes dominant above the Bragg angle, in full 
agreement with the theoretical predictions (higher bands 
become excited as the angle is further increased). 

Next, we analyze nonlinear self-action of beams with 
the inclination angles below and above the Bragg angle, 
as marked in Fig. 2. In the former case, the Bloch waves 
associated with the bottom of the first band and experi- 
encing anomalous diffraction are primarily excited. This 
results in beam self-defocusing as nonlinearity grows 10 . 
Indeed, we register broadening of the output beam with 
the increase of the laser power, see Fig. 3(a). 

Strong beam self-focusing and gap-soliton formation 
are observed when the input inclination angle is above 
the Bragg angle [see Fig. 3(b)], and the gap solitons are 
localized at narrow waveguides. Additionally, the output 
soliton position depends on the input power. Numerical 
results for the beam propagation [Fig. 4, top] agree well 
with the experimental data. In Figs. 4(a-d) we show the 
beam dynamics inside the waveguide array structure for 
various input intensities. These simulations demonstrate 
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Fig. 3. (a,b) Experimentally measured intensity distri- 
bution at the output vs. the average input power for 
incident angles marked 1 and 2 in Fig. 2, respectively. 

that the soliton velocity and, accordingly, the soliton out- 
put position are defined largely at the initial stage, where 
the counter-propagating Bloch waves of the first and sec- 
ond bands interact strongly with each other. This inter- 
band momentum exchange is predominantly responsible 
for the soliton steering, which may find applications for 
a simple realization of all-optical switching. 

In conclusion, we have generated discrete gap solitons 
in binary waveguide arrays and presented an experimen- 
tal evidence of the inter-band momentum exchange ob- 
served in the form of a power-dependent soliton steering. 
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Fig. 4. Top: Numerically calculated output field intensity 
vs. the peak input intensity (io) for the beam inclination 
20% above the Bragg angle, corresponding to Fig. 3(b). 
(a-d): Beam propagation along the binary waveguide ar- 
ray for different input intensities: (a) diffraction in linear 
regime (io = 0); (b) formation of a moving discrete gap 
soliton at io = 0.1; (c) excitation of slow or stationary 
gap solitons at io — 0.5; (d) soliton destabilization above 
a critical power threshold (Io = 1). Dashed line marks 
the boundary of experimental sample. 
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